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Abstract 

We study the ground state phase diagram and the critical properties of 
interacting Bosons in one dimension by means of a quantum Monte Carlo 
technique. The direct experimental realization is a chain of Josephson junc- 
tions. For finite-range interactions we find a novel intermediate phase which 
shows neither solid order nor superfluidity. We determine the location of this 
phase and study the critical behaviour of the various transitions. For on-site 
interaction only, we map out the phase diagram as a function of the hopping 

strength and the chemical potential. 

PACS numbers: 05.30.Jp, 67.40.Db, 67.90.+Z 
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One dimensional systems attracted strong interest over the past years. The major work 
focused on Fermionic models where a breakdown of Landau's Fermi-liquid theory opens an 
avenue for new investigations. Bosonic models deserve our attention for several reasons: 
The close analogy between Fermions and Bosons in ID, the study of enhanced quantum 
fluctuations in low dimensions, and the possibility of experimental investigation. 

A good candidate for experimental studies is 4 He, which has been investigated in the bulk 
and on surfaces 0. Another class are superconductors. The phase of the superconducting 
order-parameter is a Bosonic quantum variable. Nanofabricated circuits, especially Joseph- 
son junction networks show these quantum effects. Mott transitions have been observed 
in two dimensions || and recently also in one-dimensional systems 0. The low-energy 
Hamiltonian of a Josephson chain coincides with that of a Luttinger liquid ||. 

The interplay between the interactions and the hopping gives rise to various ground states 
which are separated by quantum phase transitions. A superconductor-insulator transition 
separates a Mott-insulating phase with localized particles from a phase with delocalized 
particles and superfluid response. The critical value of the hopping strength depends on the 
chemical potential. The inclusion of finite-range interactions enriches the picture, opening 
the possibility for various commensurate Mott-insulating phases with a solid-type ordering of 
the particles. Mean-field studies || and numerical investigations of two-dimensional systems 
showed that the transition from a Mott-insulating solid phase to the superfluid phase 
splits into two separate transitions with intermediate supersolid phases. In these supersolid 
phases superfluidity and long-range solid order coexist. No such supersolid phase has been 
found in ID 0. 

We study the quantum-phase model in ID by means of a quantum Monte Carlo method. 
For on-site interactions we determine the phase diagram with lobe-shaped Mott-insulating 
phases in good agreement with a t /[/-expansion of Freericks and Monien |ID |, as shown 



in Fig. la). With on-site and nearest-neighbour interactions the transition between solid 
phases, where the site-occupancy alternates between N and N + 1, and superconducting 
phase splits into two separate transitions as displayed by the scaling plot in Fig. 2. This 



gives rise to a novel intermediate phase which shows neither superfluidity, nor a solid order. 
We determine the location of this phase and discuss its properties. We extract the critical 
behaviour of the transitions in agreement with predictions of order-parameter descriptions 
for the onset of superfluidity fll| , |T2lj . 



The quantum-phase model describes Bosons on a lattice. Its Hamiltonian is given by 
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Two key ingredients are incorporated in this model: The interaction of the Bosons (on- 
site Uq and nearest-neighbours Ui), and the hopping t. The chemical potential fi tunes the 
number of particles on the lattice. Phase 4>i and number rij of this model are non-commuting 
operators, i.e. [rij, <pj\ = iSij. The variable i labels the lattice sites. The essential physics is 
thus dominated by the competition between phase- coherence and solid order. This gives rise 
to quantum fluctuations and quantum phase-transitions at zero temperature. The properties 
of this model are periodic in fi with a period of Uq + 2Ui. The quantum-phase model directly 
represents a chain of Josephson junctions with being the excess number of Cooper pairs 
on site i, the Coulomb interactions U, and the Josephson coupling t. A gate voltage can be 
applied to tune the chemical potential 0. The quantum-phase model is equivalent to the 
Bose-Hubbard model in the limit of a large number of Bosons per site || . 

We use the mapping of a rf-dimensional quantum model onto a (i+l-dimensional classical 
model to study the Hamiltonian of Eq. (JX|) . The ID quantum-phase model has a current-loop 
representation in 1+1 dimensions with the partition function |T3|j6|] 
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A discrete two-component current ( J x , J T ) flows on a (l+l)-dimensional space-time lattice. 
The constraint d T J T + d x J x = allows only divergence- free current loops. The imaginary 
time is discretized with the time-spacing e. The mapping is explained in detail in Ref. 0. 
The effective coupling constant K plays the role of the inverse temperature of the classical 
model and is given by K = eU /2 = In [I (et)/Ii(eb)], where J , I\ are the modified Bessel 
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functions. The time spacing e is chosen to be of the order of the inverse Josephson plasma- 
frequency, e ~ 1/y/tUo, such that the couplings in x, r-direction are isotropic. The time 
component of the current JJ can be identified as the number of excess Bosons on site i. 

The key quantities in our study are the superfluid stiffness p and the structure factor 
S(k). The former measures superfluid correlations in the system. The latter indicates 
whether the particles are arranged periodically, i.e. whether the system is in a solid ground 
state. Both quantities can be expressed in terms of the currents J on a lattice of size 
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Nearest-neighbour interactions give rise to solid phases with a finite 7r-component of the 
structure factor S w = S(k = it). 

An order-parameter description (|HJ for the onset of superfluidity predicts a dynamical 
critical exponent z = 1 in the particle-hole symmetric case at the tips of the lobes in the 
t-fi plane and a Kosterlitz-Thouless (KT) transition. For broken particle-hole symmetry 
away from these symmetry lines, z = 2 and a power law critical behaviour follows from 
these considerations. A Ginzburg-Landau description for the transition to a finite S n was 



developed in Ref. [0. In the absence of a superfluid background this action implies the 
dynamical critical exponent z = 1. 

The simulation is confined to rather small system sizes, and we use finite-size scaling for 
determining critical properties of our model. For power-law critical behaviour we use the 
scaling Ansatz |||9j 

p = L 1 -*p(L 1 ' v {t-t c ),L T /L'), (5) 

S w = L- 2 ^S(L^(t-t c ),L T /L z ), (6) 

where (3 is the critical exponent of the order parameter, v is the critical exponent for the 
correlation length, and p, S are scaling functions. 



A KT critical point can be determined using the jump of the stiffness, characteristic 
logarithmic corrections, and the behaviour of the exponentially diverging correlation length 

0- 

We also investigated the possibility of a first order transition by studying the energy 
histogram. We find no evidence for a first order transition in the phase diagram. 

We simulate the model of Eq. (0) with periodic boundary conditions using the Metropolis 
algorithm. The current loops can be divided into two classes. Local loops represent a current 
around a single plaquette on the lattice. Global loops describe a net current through the 
whole system. In each Monte Carlo sweep we try to create a local loop on each lattice site 
and global loops throughout the lattice. Each Monte Carlo run for one data point consists 
of 2 • 10 5 — 2 • 10 6 sweeps for equilibration and 10 6 — 10 7 sweeps for measurement, depending 
on the lattice size and the coupling t/U . 

We now present our numerical data, first for on-site interaction (U\ = 0). For non- integer 
values of p/U particle-hole symmetry is broken and according to Eq. (||) the scaled data 
LpQ vs. t for different lattice sizes should cross at the critical coupling t c provided L 2 /L T is 
kept constant. Fig. 16) shows the scaled superfluid stiffness vs. the coupling for p/U = 0.3. 
The curves for different lattices cross at t c jU^ = 0.207 ± 0.003. The critical exponent v is 
fitted such that plots of Lpo vs. L x l v {t — t c ) collapse onto one curve, the scaling function 
p, as shown in the inset of Fig. 16). We find v = 0.6 ±0.1 which agrees with results from 
0. At p/U = 0.4 the critical coupling is t c /U = 0.10 ± 0.005. For p/U = 0.2 we find 
t c /Uo = 0.325 ± 0.01. The critical exponent v is independent of the chemical potential 
in the range 0.2 < p/U < 0.4. At integer values of p/U , the system has particle-hole 
symmetry. Here we use the universal jump of the superfluid stiffness characteristic for a 
KT transition to determine the transition point |T^]. From the simulations with quadratic 
lattices (L, L T < 24) we obtain t c /U = 0.83 ± 0.07. 

The resulting phase diagram consists of Mott-insulating lobes with fixed integer density 
in the t-p plane, see Fig. la). The cusp- like shape of the lobes is also found in the phase 
diagram of the related Id Bose-Hubbard model analyzed in Ref. |§. In the vicinity of the 
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symmetry lines at integer values of li/Uq we find deviations from the predicted finite-size 
scaling of the superffuid stiffness of Eq. (|]). We argue that this is due to the vicinity of the 
KT transition which introduces another length scale in the scaling relation. Asymptotically 
for large systems we expect to recover a scaling according to Eq. ([5]). The solid lines in 
Fig. la) are the phase boundary obtained by a third order t/?7-expansion for the Bose- 
Hubbard model of Freericks and Monien |10| in the limit of a large number of Bosons per 



site which is in good agreement with our data. In Ref. |L0[ it is argued that the deviation 
near the tips of the lobes is due to the KT transition whose physics cannot be described in 
perturbation theory of finite order. 

When nearest-neighbour interaction is included the phase diagram is richer, including 
solid, superffuid and novel intermediate phases. We find Mott-insulating lobes where the 
site occupancy alternates periodically between N and N + 1, i.e. solid order with a finite 
structure factor S n , a superffuid phase with p 7^ 0, and a compressible intermediate phase 
with S n = and p = 0- We focus in our studies on the case of broken particle-hole 
symmetry where we can use the scaling relations of Eq. ([5]) and Eq. The dynamical 
exponent for the onset of superfluidity is thus z — 2. For the transition to a finite S n we use 
z = 1, and fit the ratio 2/3/z/ such that the scaled data cross in one point. We estimate for 
2/3 /V = 0.12 ± 0.02. Fig. 2 shows the scaled data of po and S n vs. coupling for U\/Uq = 0.4 
and fi/(Uo + 2Ui) = 0.6. We obtain for the superffuid transition t c /Uo = 0.130 ± 0.006 
and for the transition to the solid phase t c /Uo = 0.112 ± 0.005. In between there is an 
intermediate phase in which the system neither is superffuid nor solid. We address this lack 
of order at T = to the enhanced fluctuations in ID. There are corrections to the scaling 
relation of the superfluid stiffness, since the curves in Fig. 2 do not exactly cross in one 
point. A detailed investigation of these corrections, which we address again to the vicinity 
of the KT transition, is beyond our numerical resolution. 

The absence of both, superfluidity and solid order, may imply the existence of a normal 
phase of the Bosons with metallic response. A study of the superfluid stiffness as a function 
of the Matsubara frequencies and an analytic continuation can yield insight into the response 



properties [13|. Within the accuracy of our results we are not able to prove or disprove the 
existence of a finite d.c. conductivity There are general objections to a normal phase 



for Bosons at zero temperature by Leggett |T7|]. The arguments are based on the absence 
of nodes of the many body ground state wave function for Bosons in the continuum. These 
arguments are not directly applicable to our case as we study a lattice model. Furthermore, 
Leggett himself argued that the existence of a normal ground state of Bosons cannot be 
ruled out completely in ID. The question of the response of the intermediate phase remains 
open and subject to further studies. A very recent work predicts the existence of a repulsive 
Luttinger liquid intervening the solid and superfluid phase flSfl . 

In conclusion we determined the phase diagram of the quantum-phase model in ID by 
means of a quantum Monte Carlo method. For repulsive on-site interaction we mapped out 
the phase diagram as a function of the hopping strength and the chemical potential. Finite- 
range interactions give rise to new phases. The insulating phases with half-integer filling 
have solid order with a unit cell of two lattice spacings. Our simulations show an inter- 
mediate phase where both, solid and superfluid order are destroyed by the strong quantum 
fluctuations in ID. 

We would like to thank Rosario Fazio, Anne van Otterlo, Gerd Schon, and Gergely 
T. Zimanyi for stimulating discussions. This work is within the SFB195 of the "Deutsche 
Forschungsgemeinschaft" . 
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FIGURES 




FIG. 1. Results for on-site interaction, a) Phase diagram. The phase boundary separates 
the Mott-insulating (MI) phase from the superfluid phase (SF). The symbols are our Monte Carlo 



results, the solid line is the result of a third order i/C/-expansion from Ref. [10 1- b) Scaled data for 
the superfluid stiffness at fi/Uo = 0.3. The intersection of the curves gives the transition point at 
t c = 0.207 U . Array sizes L x L T : (1) 6 x 9, (2) 8 x 16, (3) 10 x 25, (4) 12 x 36. The inset shows 
the scaling of the data to a single function according to Eq. (g), with v = 0.6. 
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FIG. 2. Scaled data for the structure factor (left) and the superfluid stiffness (right) for 
Ux/Uq = 0.4 and fi/(U + 2E/i) = 0.6. In the intermediate phase for 0.112 < t/U < 0.130 the 
system shows neither solid order nor superfluidity. Array sizes L x L T for the superfluid stiffness: 
(1) 6 x 9, (2) 8 x 16, (3) 10 x 25, (4) 12 x 36. For the structure factor: L = L T = 8 (a), 10, 12, 14, 
16, 18, and 20 (b). 
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